Late-time Affleck-Dine baryogenesis after thermal inflation 
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Thermal inflation can solve serious cosmological problems such as overproduction of gravitinos 
and moduli. However, it also dilutes the preexisting baryon asymmetry. We investigate a possibility 
that Affleck-Dine mechanism works after thermal inflation and generate the baryon number at an 
acceptable level using lattice calculation. We find that a proper amount of baryon number can be 
generated for appropriate model parameters. 



PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



There are strong evidence for cosmological matter- 
antimatter asymmetry, or baryon asymmetry. Recent 
results from WMAP 3rd year data suggests [l[ 



n B 



(6.1 ±0.2) x 10 



-10 



(1) 



This value agrees well with the prediction of big-bang 
nucleosynthesis (BBN) inferred from observations of the 
light element abundances 0. Although various mecha- 
nisms are proposed in order to explain this value, it is yet 
unknown which is the right scenario. This is partly be- 
cause more or less any baryogenesis mechanism is based 
on physics beyond the standard model and partly because 
thermal history of the universe before BBN is difficult to 
probe. 

Among many candidates, perhaps the best motivated 
model as physics beyond the standard model is supersym- 
metry (SUSY) [3j. Not only it naturally solves the hier- 
archy problem, but also provides the dark matter candi- 
date. Furthermore, SUSY may play an important role in 
generating baryon asymmetry of the universe. In Affleck- 
Dine mechanism 4] the dynamics of flat directions exist- 
ing in the scalar potential of the minimal supersymmetric 
standard model (MSSM) can efficiently produce baryon 
asymmetry. 

On the other hand, the supersymmetry and supergrav- 
ity theories predicts new particles such as gravitinos and 
moduli which interact with other particles only through 
gravity and hence have long lifetimes. It is revealed that 
those particles may cause disaster for conventional cos- 
mological scenario. Among such long-lived particles, the 
cosmological problem of the gravitino has been most in- 
tensively studied by many authors. In gravity-mediation 
scenarios, gravitinos are expected to have masses of order 
100 to 1000 GeV or so, and they decay during or after 
BBN. Injected radiation and hadrons associated with the 
decay process significantly affect the light element abun- 
dances and spoil the success of BBN. In order to avoid 
this disaster, the abundance of the garvitinos should be 
small enough. Since the gravitino abundance is propor- 
tional to the reheating temperature after inflation, the 
stringent constraint on the reheating temperature is im- 
posed @, i, 0, B S M EH Q1- 

In gauge-mediation sce- 



narios, since mass of gravitino is estimated to be of order 
kcV, the gravitino becomes LSP (lightest supersymmet- 
ric particle). In this case the condition that produced 
LSP should not overdose the universe also leads to con- 
straint on the reheating temperature [l3|. Furthermore, 
in the recent articles, it was pointed out that the inflaton 
decays into gravitinos through supergravity effect, which 
results in overproduction of gravitinos even for low re- 
heating temperature [14j . 

More problematic situation arises when modulus fields 
exist. Potential of moduli is flat when SUSY is exact but 
SUSY breaking effects lift the flatness of the potential 
with curvature of order m 3 / 2 - In addition to this low en- 
ergy supersymmetry breaking, vacuum energy from in- 
flaton also breaks SUSY, which gives rise to modulus 
mass of the order of the Hubble parameter H [l5j. For 
H 3> ?7i3/2 the modulus field sits down at the poten- 
tial minimum determined by Hubble induced mass term. 
This minimum does not necessarily coincide with the true 
minimum of the potential, and naively it is expected that 
they are separated by amount of Mp(Mp: Planck mass). 
When H ~ m 3 / 2 , the modulus field starts to oscillate 
around the true minimum with typical amplitude of Mp. 
Since the coherent oscillation of modulus fields behave 
as non-relativistic matter, it soon dominates the energy 
density of the universe. As the decay rate of moduli is 
of the same order as that of gravitinos, the same cos- 
mological difficulty arises In contrast with the case 
of gravitinos lowering reheating temperature has noth- 
ing to do with moduli oscillation, so solving cosmological 
moduli problem is much more difficult. 

So far several mechanisms to avoid the cosm olog ical 
moduli problem have been proposed [l?], [H, EH l20l. Hll. 
I22I . [23l ] . Among them the most appealing solution re- 
lies on late-time entropy production to dilute the moduli. 
Here we concentrate on thermal inflation model [23|, [24| 
as a solution to the moduli problem. In the thermal 
inflation scenario a scalar field called "flaton" causes a 
mini-inflation after the primordial inflation and the sub- 
sequent decay of the flaton produce huge entropy enough 
to avoid the cosmological moduli problem (for detailed 
discussion, see 25]). It is noticed that the gravitinos are 
also significantly diluted by the thermal inflation and the 
gravitino problem is solved simultaneously. 

However, the same mechanism also dilutes the pre- 
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existing baryon asymmetry. Even Affleck-Dine mecha- 
nism can not generate the enough baryon number which 
survives from thermal inflation 26]. Thus, the baryon 
number should be regenerated after thermal inflation. 
Since thermal inflation ends at about T ~ 1 TeV, baryo- 
genesis mechanism which works at sufficiently low en- 
ergy is required. One possible candidate is electroweak 
baryogenesis [13] , but for the typical reheating tempera- 
ture after flaton decay is lower than lOOGeV, this mech- 
anism may not work. Another possibility is modified 
version of Affleck-Dine baryogenesis with LH U flat direc- 
tion which takes place after thermal inflation [28|. In 
Refs. [HI, H^] it was shown that the modified Affleck- 
Dine baryogenesis works by solving simplified dynamics 
of the scalar fields. However, actual dynamics of the 
scalar fields in this model is much complicated and it is 
not clear whether the proper amount of baryon number 
is created or not when one solves full dynamics including 
all the relevant scalar fields. In this paper, therefore, we 
study the full dynamics of the scalar fields using lattice 
simulation. We adopt the model proposed in [291 ]. 

This paper is organized as follows. In Sec[TI] we de- 
scribe our model and overview of the dynamics, in par- 
ticular how the baryon asymmetry is generated in this 
model naturally. In Sec lIIII we show the results of lattice 
simulation and constraints on the neutrino mass and fi- 
term to obtain enough baryon number. In sec lIVI we give 
our conclusions. 



II. THE MODEL 



A. fi-term and flaton field 



We start with the MSSM superpotential, 

Wmssm = y? 3 QiH u Uj +yf j Q l H d d J + yf j L i H d e j + nH u H d , 

where y's are Yukawa couplings, i is the index of genera- 
tion and SU(2) or SU(3) index is omitted. The last term 
of © is called /x-term where /i is the only dimensionful 
parameter in the MSSM. Hereafter, we also omit the gen- 
eration indices for simplicity. Furthermore, we add the 
following non-renormalizable superpotential which is re- 
sponsible for neutrino mass: 



W v = 



A„ 



2M, 



■(LH U )(LH U ), 



(3) 



Then this term explains the natural value of /i with the 
expectation value of <fi |3l| , 



V = K 



M" [ 



(5) 



where 4>q denotes the VEV of <f>. Since the required value 
of n is around electroweak scale, this relation determines 
4>q. In our scenario the most attractive choice is n = 1, 
because if we assume M ~ Mp and A p ~ 0(1) eq.© 
would give <fio ~ 10 10 GeV, which is required for flaton 
field . Therefore, in this model, the singlet field </>, 
which is introduced to solve the naturalness of /i-term, 
can be naturally regarded as flaton which leads to ther- 
mal inflation. In order to stabilize the flaton potential, 
we add the self coupling of the flaton field given by 



W A 



am' 



(6) 



where ~ 0(1). 

To summarize, our superpotential is written as 

W =y u QH u ii + y d QH d d + y e LH d e + 



2M V 



U —{LH U ){LH U 



(7) 



M 



Hi,H. 



We can forbid the other possible terms in the above su- 
perpotential using i?-parity and some discrete symmetry, 
such as Z4 symmetry, under which each field transforms 
non-trivially (29[. Furthermore, by gauging Z4 symme- 
try we can avoid the domain wall problem associated with 
discrete minima of the potential. Here we only assume 
this problem can be avoided by some mechanism which 
does not affect the dynamics we are interested in. 

Note that since the origin of neutrino mass term and 
flaton sector depend on different physics at high energy, 
cut-off scale M and M v which appear in eq.([7|) do not 
need to coincide with each other. In order to obtain a 
phenomenologically viable neutrino mass m„ < leV, The 
M v satisfies 



M v >\„x 10" GeV, 



(8) 



which suggests that M u <C Mp for the O(l) coupling. 
However, in this paper, we set M = M v ~ Mp and 
take A„ as a free parameter. Then the natural parameter 
range for \ v is expected to be rather large compared with 
X u or A<a. 



where M v denotes some cut-off scale. This term is re- 
garded as Majorana mass term for left handed neutrino 
when H u has the VEV, and also be regarded as the term 
which lifts the potential of LH U flat direction [3(|. We 
also introduce a singlet field cj> which couples with H u H d 
as 



M r - 



H,.Hd- 



(4) 



B. Scalar field dynamics after thermal inflation 

During thermal inflation all fields are trapped to the 
origin due to thermal effects, and when thermal inflation 
ends at temperature T ~ 100 GeV, the flaton field <j> and 
LH U flat direction begin to roll down to the instant min- 
imum of their potential. The full scalar potential in our 
model consists of the _F-term associated with (J7J , D-term 
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and soft SUSY breaking terms. We can parameterize the 
fiat directions LH U and H u H d as 



L 



(9) 



The F-term potential is given by Vp — J2tI> \dW/dip\ 2 , 
where tp runs all scalar fields included in the superpoten- 
tial, explicitly, 



V F =-^L{|A^ + 2X^h u h d \' z + \\JK\ 
+ \X^ 2 h d + XJ 2 h u \ 2 + |A M 0X| 2 }. 
The D-term contribution is written as 



2 |2 



V D = 9 -(\h u \ 2 -\l\ 2 -\h d \ 2 ) 2 . 



(10) 



(11) 



Finally, the soft SUSY breaking terms including soft 
masses and A-terms are 



U SB =V - ml\4>\ 2 + m 2 L \l\ 2 ~ m 2 H \h u \ 2 + m 2 Hd \h d l2 
^^ 2 hJ, 



A S X 



A„X„ 



4M 



2M 



(12) 



Note that we assume that > and m| — m H u < so 
that the flaton 4> and the flat direction LH U rolls away 
from the origin of the potential and create baryon num- 
ber after thermal inflation. Justification of these ansatz 
is rather non-trivial [32| , but it is possible to obtain valid 
parameter space [29| . In this model, as explained below, 
the most stringent constraint comes from the require- 
ment that the true vacuum must break the electroweak 
symmetry spontaneously. This implies that mass matrix 
of the Higgs fields at the origin must have the negative 
determinant, 



(H 2 -<J(H 

where B is given by 



B = A, 



2X*, 



(13) 



(14) 



from eqs. (|10[) and (|12p . Note that Hubble parameter H 
is sufficiently smaller than the typical scale appearing in 
the dynamics and negligible in the interested regime. 

Here we briefly describe outline of the dynamics of this 
model. As mentioned above, at the end of thermal in- 
flation, the flaton <fi and LH U flat direction begin to roll 
away from the origin. If the LH U direction first rolls 
away, the minimum of LH U is determined by the term 
A u X u l 2 h 2 l /2M in eq. (fl"2")) . As <j) r °Us down the potential 
and increases its field value, LH U begins to feel the posi- 
tive mass from the term \X^(j) 2 h u ^ /M 2 in eq. (fT0|) and at 
the same time, the term X^X* v (j) 2 h d l* 2 h* u / M 2 + h.c. gives 
LH U direction the angular kick. The minimum of the an- 
gular direction determined by this term is different from 



the initial one, which is source of CP violation, 
following condition 



> o 



If the 



(15) 



is satisfied, the potential for the LH U direction is stabi- 
lized at the origin and finally sit down there. However, 
it is not clear whether U(l) conserving terms dominate 
the potential at this epoch. If not, LH U direction re- 
ceives angular kick repeatedly and we lose the ability for 
predicting the resultant baryon number. Since the whole 
dynamics at this stage is quite complicated, we perform 
a numerical calculation based on lattice simulation. 



C. Reheating and baryon-to-entropy ratio 

Before going into the detailed analysis of the numeri- 
cal calculation, we comment on the reheating and entropy 
production in our model. The final reheating tempera- 
ture is determined by the flaton decay rate as 



Tr 



47r 3 g, 



(16) 



where g* denotes the effective degree of freedom at Tr. 
Let us denote the epoch at which baryon asymmetry is 
generated as ts (which is almost the same time as the 
end of thermal inflation) and at which <fi decays as t^. 
The final baryon to entropy ratio is given by 



nfl(^) n B (t B )T R 



s(t<t>) 



(17) 



Of course, it is lepton number that this process generates 
but the electroweak sphaleron effect quickly converts the 
B — L into B number [33[. 

Next we estimate the flaton decay rate. In the present 
model, 4> can decay into two Higgs bosons or higgsinos if 
kinematically allowed. The decay rate is estimated as 



■4> 



C ^T 2 



(18) 



where C is a constant of order 10 _1 ~ 10~ 2 . If we assume 
M ~ Mp, this gives Tr ~ lGeV. But note that even if 
the above decay modes are kinematically forbidden, the 
flaton can decay into two photons or gluons through the 
one loop diagrams associated with the particles in the 
thermal bath required for thermal inflation [25j. These 
processes also give decay rate that is roughly the same 
order as eq. (fT5)l . It should be noted that the natural 
range of the final reheating temperature is much lower 
than the electroweak scale. 



III. NUMERICAL ANALYSIS OF THE 
DYNAMICS 

We have investigated the full dynamics with poten- 
tial (fTU|) . (fTTj) and (fT2")l using lattice simulation. In 
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Refs. [28| and [2{|, an ad hoc large damping term was put 
by hand under assumption that the flaton field quickly 
decays through parametric resonance. In this work, we 
have not made any such artificial assumption. The only 
non-trivial input for computing the dynamics is initial 
condition. As the initial condition, we set the initial val- 
ues of all fields at around ITeV as homogeneous part, as 
is expected for thermal fluctuations at the end of ther- 
mal inflation. Our simulation is performed with one- 
dimensional lattice with 128 grid points. Initial fluctu- 
ations come from quantum fluctuation around the ho- 
mogeneous mode. We apply the method used in LAT- 
TICEEASY [3J| for including these quantum fluctuation. 
In order to eliminate unphysical effect due to large quan- 
tum fluctuations at short distance, we have cut the initial 
quantum fluctuations with mode k > m, although this 
does not much affect the result (for more detail about 
initial condition for tachyonic potential, see [35| ). 

It is found that the dynamics is rather sensitive to 
the model parameters. However, since there are many 
parameters in the model and full parameter search is be- 
yond the scope of this paper, we have performed the sim- 
ulations fixing the most of parameters as m^, — 180 GeV, 
m Hu = 700 GeV, m Hd = 800 GeV, m L = 640 GeV, 
A = 4, — 450 GeV, A v = 200 GeV, = 20 GeV, 
and arg(A A*) = -tt/4. We also take M = M P . All A- 
terms are taken to be real by field redefinition and hence 
the only remaining parameter associated with CP angle 
is arg(ApA*). The other parameters are varied in each of 
the following analysis. 

First, we show the typical motion of the field I in Fig.[T] 
when only the homogeneous mode is taken into account. 
We can see that first I rolls down to the displaced mini- 
mum and then pull back to the origin, as explained in the 
previous section. In this process, I is kicked to angular 
direction, and finally the field rotates around the origin 
with constant angular momentum (i.e., conserved lepton 
number), in the complex plane. 

From Fig. [1] it is obvious that lepton number is really 
generated. However, the present universe (or at BBN or 
recombination epoch) contains large number of causally 
disconnected regions at the era of thermal inflation. This 
means that the real baryon number is average over many 
regions with different initial values of I field. The initial 
values of I is determined by the thermal fluctuations at 
the end of thermal inflation. Thus, the initial values of 
phase of I is random and \l\ also fluctuates around ~ 
T. Taking fluctuations of I into account, we perform the 
simulations varying the initial phase of I with initial \l\ 
fixed. As for the initial value of 0, we fix its phase, but 
we have confirmed that initial angular dependence of the 
4> field does not much affect the subsequent dynamics. 
We show in Fig. [5] the relation between the initial angle 
arg(Z) and the resultant baryon asymmetry. From this 
figure it is seen that the average baryon number is really 
non-zero. Fig. [3] shows the time evolution of the baryon 
number. From this figure we can clearly see the baryon 
number is finally conserved. 
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FIG. 1: Typical motion of the field I in the complex plane. 
Here we write only zero mode. The field value is normalized 
by 10 9 GeV. We take \X„\ = 35, |A„| = 10 4 , arg(A M A;) = tt/16. 
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FIG. 2: Initial angular dependence of baryon asymmetry 
tib/s, when we take Tr — lGeV. The same parameters as 
in FigfTJare taken. 



As explained in previous section, whether or not the 
net baryon asymmetry is generated depends on CP 
phase. The phase of the flaton field takes the value that 
minimizes the potential of the angular direction deter- 
mined by the term A^X^4> 4 /M in the interested time 
scale. The relevant potential of angular direction for 
LH U comes from the term A^X u l 2 h^ initially, and sub- 
sequent dynamics depends on the terms A M A*0 2 W* 2 /i* 
and A^A* <fi 3 <fi* h^h^. Initially, the angular minimum lies 
in the direction arg(Z/i„) = 7r/2 and 37r/2. The angular 
minimum is unchanged when arg(A^A*) — it/A and 57r/4 
so the net baryon number expected to become zero. This 
is seen in Fig. [4] 

Finally, we fix arg(A M A*) = 7r/16 and choose freely the 
absolute value of A^ and X v . Each of these parameters is 
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of n ( or \fj) than about 800 GeV, although the con- 
straint (fT5|) can be satisfied, we could not get appropri- 
ate dynamics due to smallness of positive mass squared 
of LH U direction. On the other hand, the constraint (fDj)) 
invalidates /i larger than about 840GeV. As a result, only 
constrained parameter region around fi ~ 800 — 840GeV 
and 10 _3 eV< m v < 10 _1 eV can survive. 

However, we have to mention that our model has quite 
large number of parameters and the dynamics is sensi- 
tive to them. Thus the stringent constraints should not 
be taken too seriously. Complete analysis of parameter 
dependence is beyond the scope of the present paper, but 
it should be noticed that our results indicate that it is 
really possible to generate baryon asymmetry even after 
thermal inflation. 



FIG. 3: Time dependence of baryon number, t is normalized 
by (lOOGeV)" 1 . Parameters are same as FiglU and each line 
corresponds to different initial angle. Data points are reduced 
in order to make it easy to see. 
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FIG. 4: Relation between arg(A M AJ) and net baryon number. 
When CP is conserved, the net lepton number becomes zero 
as explained in the text. 



directly related to m v [eq. ©Jand \x [eq. ©]. In FigJ5]we 
plot generated baryon number (ris/s) x (Tr/GcV)" 1 in 
terms of m„ and fi. Since our choice of arg(A^A* ) = 7r/16 
corresponds to nearly maximum CP phase, the baryon- 
to-entropy ratio in Figj5] is reduced for different choice 
of arg(A M A*). Thus, the apparent large ub/s is not a 
problem. For smaller m v (or \ v ) than 10~ 3 eV, Tr 
must exceed about lOGeV even for maximum CP phase, 
which is invalid for the present model. For smaller value 



IV. CONCLUSIONS 

In this paper detailed analysis of the Affleck-Dine 
mechanism after thermal inflation has been performed. 
The present model includes only one additional gauge 
singlet which does not exist in the MSSM. This singlet 
field explains naturally both the origin of ^t-term and 
thermal inflation. This is an appealing feature when con- 
sidering the cosmological moduli problem seriously. How- 
ever, it is not trivial matter to generate enough baryon 
asymmetry consistent with late-time entropy production. 
In this paper we have demonstrated that a proper amount 
of baryon asymmetry can be generated using the lattice 
calculation for the whole dynamics without any artificial 
assumption. Although the parameter tuning is necessary 
for the present mechanism to work, It is noticed that 
since the dynamics is highly sensitive to parameters in 
the model there can be another parameter sets appropri- 
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FIG. 5: Relation among baryon number (ris/s) x 
(Tfl/GeV)" 1 , m«(eV) and /it(GeV). Region below the hori- 
zontal line requires Tr > lOGeV to obtain desired value of 
n B /s. 
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ate for our purpose which we could not find. 

We stress that if the late time entropy production takes 
place and dilute the preexisting baryon number density, 
the modified Affleck-Dine process considered in this pa- 
per is the only known mechanism to re-create the baryon 
number at very low energy scale after late-time entropy 



production, ft is interesting that this mechanism can 
work with quite low reheating temperature even below 
the electroweak scale. This opens up a new possibility 
for such a low-energy baryogenesis, and may be tested in 
future accelerator experiments. 
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